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We consider the fractional generalizations of the phase volume, volume element and Poisson 
brackets. These generalizations lead us to the fractional analog of the phase space. We consider 
systems on this fractional phase space and fractional analogs of the Hamilton equations. The 
fractional generalization of the average value is suggested. The fractional analogs of the Bogoliubov 
hierarchy equations are derived from the fractional Liouville equation. We define the fractional 
reduced distribution functions. The fractional analog of the Vlasov equation and the Debye radius 
are considered. 
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I. INTRODUCTION 

Derivatives and integrals of fractional order have found 
many applications in recent studies of scaling phenom- 
ena [1 U H II H In Ref. 0, coordinate fractional 
derivatives in the Fokker-Planck equation were used. It is 
known that Fokker-Planck equation can be derived from 
the Liouville equation ||. The Liouville equation is de- 
rived from the normalization condition and the Hamilton 
equations 0. In the Hamilton equations dqk/dt — Pk/m, 
dpk/dt = Fk(q,p) we have only the time derivatives. The 
usual normalization condition leads to the usual (non- 
fractional) Liouville equation. If we would like to de- 
rive the fractional Liouville equation then we must use 
a fractional normalization condition. In Ref. |lfj|. the 
fractional Liouville equation was derived from fractional 
normalization condition. 

The natural question arises: What could be the phys- 
ical meaning of the fractional normalization condition? 
This physical meaning can be the following: the system 
is to be found somewhere in the fractional phase space. 
The fractional normalization condition can be considered 
as a normalization condition for the distribution function 
in a fractional phase space. In order to use this interpre- 
tation we must define a fractional phase space. The first 
interpretation of the fractional phase space is connected 
with fractional dimension space. The fractional dimen- 
sion interpretation follows from the formulas for dimen- 
sional regularizations . This interpretation was suggested 
m Ref. ha. In this paper we consider the second inter- 
pretation of the fractional phase space. This interpreta- 
tion follows from the fractional measure of phase space 
[Io| that is used in the fractional integrals. The frac- 
tional phase space is considered as a phase space that is 
described by the fractional powers of coordinates and mo- 
menta (q%,p%). In this case, the fractional normalization 
condition for the distribution function and the fractional 
average values are considered as a condition and values 
for the fractional space. In general, these systems are 
non-Hamiltonian dissipative systems for the usual phase 



space (q,p). 

It is known that Bogoliubov equations can be derived 
from the Liouville equation and the definition of the av- 
erage value 0111 m 0. In Ref. [ljj, the fractional 
Liouville equation is derived from the fractional normal- 
ization condition. In this paper we define the fractional 
analog of the average value and reduced distributions. 
Then we derive the fractional Bogoliubov equations from 
the fractional Liouville equation and the definition of the 
fractional average value. We derive the fractional analog 
of the Vlasov equation and the Debye radius. 

In Ref. OH, the Ricmann-Liouville definition of the 
fractional integration and differentiation is used. There- 
fore in this paper we use this definition of fractional in- 
tegration and differentiation. 

In Sec. II, we define the fractional phase space vol- 
ume. The fractional phase volume element for the frac- 
tional phase space is considered. We define the fractional 
analog of the Poisson bracket. In Sec. Ill, we consider 
the fractional systems. We discuss the free motion of 
the fractional systems, the fractional harmonic oscilla- 
tor and fractional analog of the Hamiltonian systems. In 
Sec. IV, the fractional average values and some notations 
are considered. We define the reduced one-particle and 
two-particle distribution functions. In Sec. V, the frac- 
tional Liouville equation for n-particle system is written. 
We derive the first fractional Bogoliubov equation from 
the fractional average value and the fractional Liouville 
equation. The second fractional Bogoliubov equation is 
considered. In Sec. VI, we derive the fractional analog of 
the Vlasov equation and the Debye radius for fractional 
systems. Finally, a short conclusion is given in Sec. VII. 



II. FRACTIONAL PHASE VOLUME AND 
POISSON BRACKETS 

A. Fractional Phase Volume for Configuration 
Space 

Let us consider the phase volume for the region such 
that x G [a; b\. The usual phase volume is defined by 
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Vi= dx = 



dx 



dx. 



(1) 
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The fractional integrations are defined ^(| by 

n- 1 r dx 

a+ T(a)J a {y-xY-^ 



4-1 



f r dx 
T(a)J y (x-yY-*- 



Using these notations, we get Eq. in the form 

V 1 = 7 a 1 + l + / b 1 _l. (2) 

The fractional generalization of this phase volume can be 
defined by 



v a = i« + i + i£_i 



(3) 



The fractional phase volume integral can be written 
in the form 



1 + A° 



a pb-y 



r(«) Jo 



(4) 



where A is defined by 



V -a 
b-y' 



In order to have the symmetric limits of the phase vol- 
ume integral we can use the following equation: 

V n = 1 t > I dfi a (x), (5) 



where 



djjL a (x) 



ib-y) 



l dx dx° 



F(a) aT(a) 



(6) 



Here we use the following notation for fractional power 
of coordinates: 



x a = f3{x)(x) a = sgn(x)\x\ a , 



(7) 



where f3(x) — [sgn(x)] a . The function sgn(x) is equal 
to +1 for x > 0, and this function is equal — f for x < 0. 



B. Fractional Phase Volume for Phase Space 

Using Eq. JSJ, we have the phase volume for the two- 
dimensional phase space in the form 



Va 



where 



(i + A?) (I + A") /■+(»-») r+(p»-v) dq a A dp 



2 2 J-{<v,-y) J-(p b -v') ( aT ( a )) 



(iq Q A dp" = e^lgpl - 1 ^ A dp. 



(8) 
(9) 



The fractional measure for the region B of 2n- 
dimcnsional phase space can be defined by the equation 

/j, a (B) = V a = / g{a)dn a {q,p), (10) 

JB 

where dfj, a (q,p) is a phase volume element, 

(oT(a)) 2 

and (7(a) is a numerical multiplier, 



(11) 



1 ™ 



fc=i 



If the domain B of the phase space is defined by qu 6 R 1 
and pfc S R , then gk(ct) = 1. If this domain is defined 
by g/c € [<7 a /c;<7&fc] andp fc e [p a /=;Pf>fe]i then 



flfc (a) = 1 + 



qbk - Vk 

Vk - Qak 



1 + fp^zAy\ (12) 



It is easy to see that the fractional measure depends on 
the fractional powers (q a ,p a ). 



C. Fractional Exterior Derivatives 
In Eq. IjSJl, we use the usual exterior derivative 

2» 



d = dxh 



k=l 



d 

dx k ' 



(13) 



Obviously, this derivative can be represented in the form 

2 n 



d = Ydx%^- 



k=l 



dx° 



(14) 



where x a is defined by Eq. J7J. 

Note that the volume element of fractional pha se sp ace 
can be realized by fractional exterior derivatives |l5|. 

d ° = S dqt - »))- + S (d(Pk-y' k )r ' 



in the following form: 
4 



dq a Adp a = 



1 N- 1 



r 2 (2-a) T 2 (l - a) 



) {qp) a - 1 d a qAd a p. 



D. Fractional Poisson Brackets 

We can define the fractional generalization of the Pois- 
son brackets in the form 

d a A d a B d a A d a B \ 
i ' B|(Q) ~ ^ W) Q (d Pk ) a ~ (dp k ) a {dq k ) a r 

(15) 
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where we use the fractional derivatives ^(|. It is known 
fl(| that the derivative of a constant need not be zero, 



d a l 



{dx) a T(l - a) 



(16) 



This equation leads us to the correlation between coor- 
dinates q and momenta p in the form 



d a p 



pq 



d a q 



qp 



(dq) a r(l-a)' {dp) a r(l-a)' 



(17) 



It is easy to see that q and p are not independent variables 
in the usual sense. Therefore the fractional analog i|15|) 
of the Poisson brackets is not convenient. 

In the general case (n > 2), the Poisson brackets 
{qk,qi}( a ) an d {pk,Pi}( a ) (with k ^ I) are not equal to 
zero. If we consider n = 2, then we have 

{<7i,92}(a) = -j2(a)qiq 2 [(qiPiy a - (gs^a) - "], (18) 
where the coefficient 72(a) is defined by the equation 
1 2 



72(a) = 



r 2 (l-a) r(l - a)r(2 - a) ' 



(19) 



Moreover, the Poisson brackets with the unit are not 
equal to zero. Using Eq. JTBJl, we get 



{l,q} {a) = l2 (a)q 1 - a p a . / 



{1,p}(q) = -72(a)? a p 1 



(20) 



(21) 



Therefore the fractional Poisson brackets (|15|) are not 
convenient. 

We can use the fractional power Poisson brackets: 
where we use notations {7J). Obviously, the relations 

{<?fc,<?f}a = 0, {pg,p?}a=0, {g?,P?} a = <y«, 

are realized for these Poisson brackets. Therefore the 
fractional Poisson brackets are more convenient. 

E. Fractal Dimension of Space 

The interpretation of the fractional phase space is con- 
nected with the fractional measure of phase space. The 
parameter a defines the space with the fractional phase 
measure (|1(JH and l|ll|) . It is easy to prove that the veloc- 
ity of the fractional phase volume change is defined by 
the equation 

d\^ f 

= J Cl a (q,P,t)g{a)dfi a (q,p), 



where 



r dq? a i 1 r a ^P? 1 

{-JT^Pt } a + {q t ,-jr}c 



dt 



dt 



(23) 



Equation. I|23|) is proved in Ref. . The form of the 
omega function allows us to consider a different class of 
the systems that are described by the fractional powers 
of coordinates and momenta. 

The interpretation of the fractional phase space can be 
derived from the fractional measure that is used in the 
fractional integrals. The interpretation of the fractional 
phase space can be connected with fractional dimension. 
We have two arguments for this point of view. 

1. Let us use the well-known formulas for dimensional 
regularizations |l7j : 



p{x)d n x 



2W 2 
r(n/2) Jo 



p{x)x n ' x dx. (24) 



Using Eq. 1)24(1 , we get that the fractional integral can be 
considered 10j as an integral for the fractional dimension 
space 



r(q/2) 

2W 2 r(a) 



p(x, t)d a x 



(25) 



up to the numerical factor r(a/2)/(27r Q / 2 r(a)). 

2. Let us consider the well-known definition of the 
fractal mass dimension. The equations that define the 
fractal dimensions have the passage to the limit. This 
passage makes difficult the practical application to the 
real fractal media. The other dimensions, which can be 
easily calculated from the experimental data, are used 
in the empirical investigations. For example, the mass 
fractal dimension 0, 0] can be easy measured. 

The properties of the fractal media like mass obeys a 
power law relation, 



M(r) = kr 



(26) 



where M is the mass of fractal medium, r is a box size 
(or a sphere radius), and D m is a mass fractal dimension. 
The amount of mass of a medium inside a box of size r 
has a power law relation (|2"B*|l . 

The fractal dimension can be calculated by box count- 
ing method which means drawing a box of size r and 
counting the mass inside. To calculate the mass fractal 
dimension, take the logarithm of both sides of Eq. 12611 : 

ln(M) = D m ln(r) + Ink. 

The log- log plot of M and r gives us the slope D m , the 
fractal dimension. When we graph ln(M) as a function 
of ln(r) , we get a value of about D m which is the fractal 
dimension of fractal media. 

The power law relation (|26|l can be naturally derived 
by using the fractional integral. In order to describe the 
fractal media, we suggest to use the space with fractional 
measure. 
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Let us consider the line distribution of the mass. If 
we consider the mass of the homogeneous distribution 
(p = const) in the ball region W with radius r, then we 
have 

/+r rr 
p{x)dx = 2p dx = 2r\ (27) 
-r JO 

In this case, D m = 1. Let us consider line mass distribu- 
tion in the fractional space. In that case, a ball of radius 
r covers a mass 



M a {r) 



2pr a 
aT(a) 



This equation can be proved by the fractional generaliza- 
tion of Eq. H27JI in the form 

M a (r) = f piO^- 1 ^ = P V a {r) 2/ "" 



r (o0 Jo 



aT(a) 



(28) 

The initial points in the fractional integral (@J are set 
to zero, and a = —r, b = r. The fractal dimension of 
particle system and fractal media is defined as the expo- 
nent of r in the growth law for mass M(r) or number of 
particles, 



n(r) = M(r)/m = (k/m)r l 



(29) 



Here m is a particle mass. Thus we see that the fractal 
dimension D m of particle system (in the fractional space) 
is a, i.e., D m — a. Therefore the space with fractional 
measure Q can be considered as a space with fractal 
dimension D m — a. 

As the result the space with fractional measure iQJ) can 
be used to describe the particle systems and medium with 
non-integer mass dimension. 



III. FRACTIONAL SYSTEMS 
A. Equations of Motion 

In Sec. II we prove that the fractional measure 
depends on the fractional powers (q a ,p a ), and Poisson 
brackets (|22f) with fractional powers are more convenient. 
Therefore we can consider a different class of mechanical 
systems that are described by the fractional powers of co- 
ordinates and momenta. We can consider the fractional 
power of the coordinates as a convenient way to describe 
systems in the fractional dimension space. 

Let us consider a classical system with the mass M. 
Suppose this system is described by the dimensional co- 
ordinates q k and momenta p k that satisfy the following 
equations of motion: 



dq k p k dq k . , _ 



dt 



M 



dt 



(30) 



Let qo be the characteristic scale in the configuration 
space; po be the characteristic momentum, Fq be the 



characteristic value of the force, and to be the typical 
time. Let us introduce the dimensionless variables 

q~k Pk , i f k 

qk = — , Pk = —, t=—, F k = — . 
qo Po t F a 

Here and later we use q k and p k as dimensionless vari- 
ables. Using Eq. (|30|l for dimensional physical variables 
(<7, p), we get the equations for dimensionless variables, 

dq k p k dq k 

-rr = — , ~ff = AF k (q,p,t), (31) 
at m dt 

where m = Mqo/topo is the dimensionless mass, and 

toF 



A 



Po 



(32) 



Using the dimensionless variables (q,p,t), we can con- 
sider the fractional generalization of Eq. I|31|) in the form 



dq% _ Pt dqt 



AF k (q a ,p a ,t), 



dt m dt 
where we use the following notations: 

qt=m(lk) a =sgn(qk)\qk\ a , 
p a k =[}(p){PkT = sgn{p k )\p k \ a . 



(33) 



(34) 



(35) 



Here k = 1, ...,n, and j3(x) is defined by Eq. J7J. 

A system is called a fractional system if the phase space 
of the system can be described by the fractional powers 
of coordinates (|34|l and momenta (|35|) . 

The fractional phase space can be considered as a phase 
space for the fractional systems. This interpretation fol- 
lows from the fractional measure that is used in the frac- 
tional integrals. 

We can consider the fractional systems in the usual 
phase space (q,p) and in the fractional phase space 
(q a ,p a ). In the second case, the equations of motion for 
the fractional systems have more simple form. Therefore 
we use the fractional phase space. The fractional space 
is considered as a space with the fractional measure that 
is used in the fractional integrals. 

The fractional generalization of the conservative 
Hamiltonian system is described by the equation 



dt ~ ^ ~ 1 j 



dH 

dpi 



dp% 
dt 



dH 

dq%' 



where H is a fractional analog of the Hamiltonian. Note 
that the function H is the invariant of the motion. Using 
the fractional Poisson brackets (|22|l . we have 



d € _ r n a m dpi 



{Pk,H} a . 



(37) 



Here we use Poisson brackets (|22H . These equations de- 
scribe the system in the fractional phase space (q a ,p a ). 
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For the usual phase space {q,p), the fractional Hamilto- 
nian systems are described by the equations 



dgk 
dt 



(qkPk) 1 



OH 



dpk 



dpk 
dt 



{ikPk) 



dH 



dq k 



(38) 



The fractional Hamiltonian systems are non-Hamiltonian 
systems in the usual phase space (q,p). A classical system 
(in the usual phase space) is called Hamiltonian if the 
right-hand sides of the equations 



dqk f s dp k . 
-rr = 9k{q,P), -TT = fk{q,P) 

satisfy the following Helmholtz conditions |2C 



(39) 



%fc__% =0 ^9k_^Il^Q ^J±_^A =Q (40) 
dpi dpk ' dqi dpk ' dqi dqk 

It is easy to prove that the Helmholtz conditions are not 
satisfied. Therefore the fractional Hamiltonian system 
(|38|l is a non-Hamiltonian system in the usual phase space 
(q,p). The fractional phase space allows us to write the 
equations of motion for the fractional systems in the sim- 
ple form 

If we have dq^/dt = pt/m, then the fractional analog 
of the Hamiltonian can be considered in the form 



a+f) 
a Pk 

m(a + (3) 



U(q). 



(41) 



If j3 — a, then we have the fractional analog of the 
kinetic energy T = p 2a /2m. 

The omega function for the system H39JI in the usual 
phase space (q,p) is defined by the equation 



CI 



fc=i 



dpk 



(42) 



If the omega function is negative Cl < 0, then the sys- 
tem is called a dissipative system. If fi ^ 0, then the 
system is a generalized dissipative system. For the frac- 
tional Hamiltonian systems (|38|) , the omega function i|42|) 
is not equal to zero. Therefore the fractional Hamiltonian 
systems are the general dissipative systems in the usual 
phase space. 

It is known that the non-Hamiltonian and dissipative 
systems (for example, F = —jp) are not invariant under 
the Galilean transformation. In general, the fractional 
systems are non-Hamiltonian systems for the usual phase 
space (q,p). Therefore the Galilean transformations for 
the equations of motion are not considered in this paper. 
The fractional analogs of the Hamiltonian systems can 
be invariant under the fractional generalization of the 
Galilean transformation. 

It is not hard to prove that the fractional systems i|41|) 
are connected with the non-Gaussian statistics. Classi- 
cal dissipative and non-Hamiltonian systems can have 
the canonical Gibbs distribution as a solution of the 



stationary Liouville equations [2l|. Using the methods 
|21j. it is easy to prove that some of fractional dissipa- 
tive systems can have fractional Gibbs distribution (non- 
Gaussian statistic) 



p(q,p) = exp [F - H(q,p)]/kT, 



(43) 



as a solution of the fractional Liouville equations [To| . 
The interest in and relevance of fractional kinetic equa- 
tions is a natural consequence of the realization of the 
importance of non-Gaussian statistics of many dynam- 
ical systems. There is already a substantial literature 
studying such equations in one or more space dimensions. 

B. Free Motion of Fractional System 

Let us consider the free motion of the fractional system 
that is defined by the following equations: 



dq a p a dp a 



dt 



dt 



= 0. 



(44) 



The solutions of these equations with the conditions 
q(0) = q\ and p(0) — p\ have the form 



q a (t) =m' 1 p^t + q*, p(t) = Pl . 



(45) 



We can conclude that the free motion of the fractional 
system is described by 



g(*)~(t-ti) 1/c 



(46) 



where the parameter t\ is defined by q\. 

Let us consider the special cases of the parameter a: 
0.5, 1/3, 0.6. 

1. If the parameter a is equal to 0.5, then we have the 
solution of Eq. 144|) in the form 

q(t) = z( qi , Pl ,t){at 2 + bt + c), (47) 

where we use the following notations 

z{qi,Pi,t) = sgn{y/\jn\t - my/\q~[\), (48) 

a = \px\/m, b = 2^\q lPl \/m, c=|gi|. (49) 

2. If the parameter a is equal to 1/3, then we have the 
solution of Eq. I|44|l in the form 

o 2/3 1/3 o 1/3 2/3 

q(t) = ^ + 3 J^^t 2 + 3 -^^t + qi . (50) 

3. If the parameter a is equal to 0.6, then we have the 
solution of Eq. <|44|l in the form 



5/3 



(51) 



If ft = 0, then q(t) - t 5 / 3 . 
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For the usual phase space (g, p) , the equations of mo- 
tion for the fractional system can be represented in the 
following form 1)39(1 . For the free fractional system, we 
have 



dgk _ 

dt am 



q l-a p a 



dpk 
dt 



= 0. 



(52) 



The omega function for the free fractional system is equal 
to the following function: 



1 



a 



am 



(53) 



In general, this function does not equal to zero and the 
phase volume of the usual phase space changes. Using 
Eq. we get 



^ 1 — a , ,i 



(54) 



where t\ — mq" /p". If q\ — 0, then the omega function is 
proportional to 1/t. Therefore the velocity of elementary 
phase volume change for the free motion in the usual 
phase space (q,p) is inversely proportional to the time. 
We suppose that the initial momentum is not equal to 
zero pi 7^ 0. 

For the fractional phase space (q a ,p a ), we define 01 
the omega function Q a) in Eq. 1(23(1 . This "fractional" 
omega function is equal to zero for the free motion of 
fractional system. Therefore the fractional phase space 
(q a ,p a ) is more convenient than usual phase space (q,p). 



C. Fractional Harmonic Oscillator 

Let us consider the fractional harmonic oscillator, 
which is defined in Ref. by the equations 



dq a p a dp a 9 „ 

-jr = — i -jt = -rnwV, 
at m at 



(55) 



where w is a dimensionless variable. The solutions of 
these equations of motion have the form 

q a (t) = a cos(ujt + if), p a (t) = —mioa sin(ujt + if), 

(56) 

where parameters a, if are defined by 



Pi 



uj^m 

If a = 0.5, then wc have 



mujq^ 



(57) 



q(t) = sgn(cos(u>t + fi)j a 2 cos 2 {uit + if). (58) 

If a = 1/3, then we have q(t) = a 3 cos 3 (ujt + if). The 
fractional harmonic oscillator has the following integral 
of motion: 



p 2a mui 2 q 2a 

H = — 1 = const. 

2m 2 



(59) 



This function can be considered as a fractional analog of 
Hamiltonian for Eq. 1)36(1 . 

For the fractional harmonic oscillator the function f2 Q 
is equal to zero in the fractional phase space (q a ,p a )- 
Therefore this system is a conservative Hamiltonian non- 
dissipative system in the fractional space. If we use 
the usual phase space to describe the fractional har- 
monic oscillator, then this system is conservative non- 
Hamiltonian dissipative system. Note that the conser- 
vative non-Hamiltonian systems are considered in Ref. 

E2 

In the usual phase space (<?, p) the equations of motion 
for fractional oscillator have the form 



dq 
~dt 



1 



dp 
~dt 



-q a p 1 ~ a . 



(60) 



The omega function for the usual phase space (q,p) is 
defined by Eq. l(l2*|l in the form 



n 



p (p — m ui q ). 



am 



(61) 



The elementary phase volume of the usual phase space 
changes. The velocity of this change is equal to the omega 
function. Substituting the solution l(56|) into Eq. I|61|) , we 
have 



CI = 



1 



a 



a 



2u) cot 2{ujt + if), 



(62) 



where we use sin 2 j3 — cos 2 '(3 — —cos2f3 and cot(3 — 
cos(3 1 'sin/3. Therefore the fractional harmonic oscillator 
is a general dissipative system in the usual phase space 

For the fractional phase space (q a ,p a ), the fractional 
harmonic oscillator is conservative nondissipative Hamil- 
tonian system. Therefore the fractional phase space 
(q a ,p a ) is more convenient than usual phase space (q,p). 

The question arises: What are the fundamentals (dif- 
ferent from Hamilton principle), which can lead to the 
system of dynamic equations ((60|) ? For the fractal me- 
dia the harmonic oscillator is defined by Eq. ((55(1 . For 
the usual phase space (q, p) , this equation has form 16011 . 
Therefore the system ((60(1 can be considered as a system 



dq 
~dt 



P_ 

M ' 



dp 

a 



= -M{Lo/a) 2 q, 



(63) 



where the dimensionless mass M(q,p) = ma^p/q) 1 a 
satisfies the scaling relation 

M(X iq , X 2P ) = (A 2 /A 1 ) 1 - Q M( g ,p). 

This property can be formulated in the following form. If 
we consider the scale transformation of the characteristic 
values in the form 

qo -> goAi, Po -> P0/A2, 

then we have transformation of the mass 

M -» (A 2 /Ai) 1 ~ Q M. 
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In the general case, this scaling law can be described by 
the renormalization group approach 23] . 

Note that the system (|fc>3l) is non-Hamiltonian system. 
We consider the fractional phase space form l|55ll of Eq. 
(|60|l as a more fundamental. The fractional harmonic os- 
cillator is an oscillator in the fractional phase space that 
can be considered as a fractal medium. Therefore the 
fractional oscillator can be interpreted as an elementary 
excitation of some fractal medium with noninteger mass 
dimension. 



IV. FRACTIONAL AVERAGE VALUES AND 
REDUCED DISTRIBUTIONS 

A. Fractional Average Values for Configuration 
Space 

Let us derive the fractional generalization of the equa- 
tion that defines the average value of the classical observ- 
able A(q,p). 

The usual average value for the configuration space 



< A >i= / A(x)p{x)dx 



(65) 



D. Curved Phase Space 

The fractional system with the fractional analog l|41|) of 
the Hamiltonian can be considered as a nonlinear system 
with 



n 1 

H(q a ,p a ) = J2 ^9ki{q,p)PkPi + U(q). (64) 
fc=l 



Note that this fractional Hamiltonian (16411 defines a non- 
linear one-dimensional sigma-model [24j,[25j in the curved 
phase space with metric gki{q,p) = m^P^^^dki- This 
means that we use the curved phase space. Note that this 
fractional Hamiltonian is used in equations of motion 138|) 
that define the non-Hamiltonian flow in the usual phase 
space. 

The curved phase space is used in the Tuckerman 
approach to the non-Hamiltonian statistical mechanics 
lilllllilEillillSElElSlli- In their ap- 
proach the suggested invariant phase space measure of 
non-Hamiltonian systems is connected with the metric 
of the curved phase space. This metric defines the gener- 
alization of the Poisson brackets. The generalized (non- 
Hamiltonian) bracket is suggested in Refs. |3ll l32j|. For 
these brackets, the Jacobi relations will not be satis- 
fied. This requires the application of non-Lie algebras 
(in which the Jacobi identity does not hold) and ana- 
lytic quasi-groups (which are nonassociative generaliza- 
tions of groups). In the paper j3f|, we show that the ana- 
logues of Lie algebras and groups for non-Hamiltonian 
systems are Valya algebras (anticommutative algebras 
whose commutants are Lie subalgebras) and analytic 
commutant-associative loops (whose commutants are as- 
sociative subloops (groups)). Unfortunately, non-Lie al- 
gebras and its representations have not been thoroughly 
studied. Nevertheless the Ricmannian treatment of the 
phase space is very interesting. This approach allows us 
to consider the connection between the fractal dimen- 
sional phase space and non-Lie algebras of the vector 
fields in the space. 



can be written in the form 

/y roo 
A(x)p(x)dx + / A(x)p(x)dx. (66) 
-oo Jy 



Using the notations 



mm = 



1 f y f(x)dx 



r(«) 7_oo (z/-*) 1 -"' 



(i-f)(y) = 



T(a) Jy (x - y) 



f(x)dx 

l-a ' 



average value (|66|l can be rewritten in the form 

< A > 1= (Il Ap)(y) + (lUp)(y). 

The fractional generalization of this equation is defined 
by 

< A > a = (I$Ap)(y) + {I a _Ap){y). (67) 
We can rewrite Eq. I|67|) in the form 



< A > a = - / ((Ap)(y -x) + {Ap)(y + x))dp a (x), 



where we use 



dp a (x) 



1 dx dx a 



T(a) aT(a) 
and x a is defined by Eq. J7J). 



(68) 
(69) 



B. Fractional Average Values for Phase Space 

Let us introduce some notations to consider the frac- 
tional average value for phase space. 

1. Let us define the tilde operators T Xk and T[k]. The 
operator T Xk is defined by 

T Xk f(..-,x k , ...) = -(f(...,x' k -x k , ...)+f(...,x' k +x k , ...))■ 



8 



This operator allows us to rewrite the functions 
- (Atf -q,p'~ p, t)p(q' -q,p'- p, t) + 

+A(q' +q,p'- p, t)p(q' + q,p' - p, t) + 
+A(q' -q,p' + p, t)p(q' - q,p' + p, t) + 

+A(q +q,p + p, t)p(q +q,p' + p, t)j 

in the simple form 

T q T p (A(q,p,t)p(q,p,t)). 

Let us consider k particle that is described by general- 
ized coordinates qks and generalized momenta pks , where 
s = 1, ...,m. The operator T[k] is defined by the relation 
T[k] = Tq kl Tp kl ...T qkrr T Pkm . For the n-particle system 
phase space, we use the operator T[l, ...,n] — T[l]...T[n]. 

2. Let us define the integral operators Z? and j a [fc]. 
The operator I" k is defined by the equation 



Ix k f{ x k) = / f{x k )dp a (x k ). 



(70) 



In this case, fractional integral (|68fl . which defines the 
average value, can be rewritten in the form 

< A > a = I«T x A{x)p(x). 

Let us define the phase space integral operator I a [k] 
for k particle by I a [k] = Jq k Jp kl -Jq km ip km , i-e., we use 

I a [k\f (<lk,Pk) = / /(qfc,Pfc)eW(qfc,P*0- (71) 



Here dp a ((\ k , p k ) is an elementary 2m-dimensional phase 
volume that is defined by the equation 

dp a (q k , Pk ) = (aT(a))- 2m dq k * 1 Adp k * 1 A...Adq% m Adp% m . 

For the n-particle system phase space, we use the integral 
operator I a [l, ...,n] = i a [l]..J a [n]. 

3. Let us define the fractional analog of the average 
values < A > for the phase space for n-particle system. 
Using the suggested notations, we can define the frac- 
tional average value by the relation 



< A> a =i a [l,...,n]T[l,...,n]A Pn . 



(72) 



In the simple case (n = m = 1), the fractional average 
value is defined by the equation 



< A > Q = 



OO rOC 



dp, a (q,p) TqT p A(q,p)p(q,p). (73) 



OO J — OO 



Note that the fractional normalization condition 01 
is a special case of this definition of the average value 



C. Reduced Distribution Functions 

In order to derive a fractional analog of the Bogoliubov 
hierarchy equations we must define the reduced distribu- 
tions. 

Let us consider a classical system with fixed num- 
ber n of identical particles. Suppose k particle is de- 
scribed by the dimensionless generalized coordinates q ks 
and generalized momenta p ks , where s = l,...,m. We 
use the following notations qfc = (q k i, q km ) and pfe = 

(pfcl, -,Pkm)- 

The state of this system can be described by dimen- 
sionless n-particle distribution function p n in the 2mn- 
dimensional phase space 

p„(q,p,<) = p(qi,pi, ...,q„,p„,i). 

We assume that distribution function is invariant un- 
der the permutations of identical particles [3(| : 

p(...,q k ,p k , ...,qi,p;, ...,i) = p(...,qi,pi, ...,q fc ,p fc , 

In this case, the average values for the classical observ- 
ables can be simplified. We use the tilde distribution 
functions 

ft,(q,p,t) =T[l,„.,n]p n (q 1 p,t), (74) 

and the function p\ that is defined by the equation 

pi(q,p,t) = p(qi,pi,i) = i a [2,...,n]p n (q,,p,t). (75) 

This function is called one-particle reduced distribution 
function. The function is defined for the 2m-dimensional 
phase space. Obviously, that one-particle distribution 
function satisfies the normalization condition llOl 



J Q [l]p 1 (q,p,t) = l. 



(76) 



< 1 > a = 1. 



Two-particle reduced distribution function p>2 is de- 
fined by the fractional integration of the n-particle dis- 
tribution function over all q k and p kl except k = 1,2: 

P2(q,P,f) = p(qi,Pi,q 2 ,P2,i) = I a [3, n]p„(q, p, t). 

(77) 



V. FRACTIONAL LIOUVILLE AND 
BOGOLIUBOV EQUATIONS 

A. Fractional Liouville Equation 

The fractional generalization of the Liouville equation 
is derived in Ref. ^jj. Let us consider the Hamilton 
equations for n-particle system in the form 

" C ~G k s {q a ,P a ), ^ = AF*(q a ,p a ,t). (78) 



dt 



dt 
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The evolution of n-particle distribution function p n is 
described by the Liouville equation. The fractional Li- 
ouville equation 0] for n-particle distribution function 
has the form 



dpn 

dt 



&aPn = 0. 



(79) 



This equation can be derived |To) from the fractional nor- 
malization condition 

I a [l,...,n]p n (q,p,t) = l. (80) 

In the Liouville equation d/dt is a total time derivative 
n,m , „ n,m , „ 

A ~ dt dt dq ks ^ (ft op fcs ' 

Using Eq. J7SJ), this derivative can be written for the 
fractional powers in the form 

The a-omega function is defined by the equation 

71,771 

= E (i G tPk s U+A{qt,F^} a ). (82) 

fc,a=l 

Here we use the following notations for the fractional 
Poisson brackets: 



/ a m _ ( dA dB 



5^ 95 as 



(83) 



Using 17911 , (|82(l and (|81|l , we get the Liouville equation 
in the form 



at 



(84) 



where A„ is Liouville operator that is defined by the equa- 
tion 



a - ( d{G h sPn) , A d{F*p n ) \ 



fe,s=l 



OP? 



B. First Fractional Bogoliubov Equation 

The Bogoliubov equations [HEEE! 

are equations 

for the reduced distribution functions. These equations 
can be derived from the Liouville equation. Let us derive 
the first fractional Bogoliubov equation (|93|) from the 
fractional Liouville equation Q84[l. 

In order to derive the equation for the function p\ we 
differentiate Eq. (J7SJ, which defines one-particle reduced 
distribution: 

f = ^ P ,...,»]P„ = /-P,.,n]f. 



Using the Liouville equation (|84|l for n-particle distribu- 
tion function p n , we have 



dpi 
dt 



7 Q [2,...,n]A„p„(q,p,i). (86) 
Substituting Eq. in Eq. (JSHJ), we get 

dpi fQro ^ (d(G k s ~ Pn ) l A d{FipnY 



fc,s=l 



9pL 



(87) 



Let us consider in Eq. (|87|l the integration over q ks 
and pfc s for k-particle term. Since the coordinates and 
momenta are independent variables, we can derive 



d 



l 



Here we use the condition 



G k s p n 



= 0. 



lim p n = 0, 

q ks —t±oo 



(89) 



which follows from the normalization condition. If limit 
(|59"|) is not equal to zero, then the integration over phase 
space is equal to infinity. Similarly, we have 

Then all terms in Eq. (|87|l with k = 2, n are equal to 
zero. We have only term for k — 1. Therefore Eq. I|87|l 
has the form 



dpi 
dt 



5>[ 2 ,...,»i(ffif 



A«3hl). « 



Since the variable qi is an independent of q2, ...q n and 
P2, ...p n , the first term in Eq. I|90|) can be written in the 
form 



d{G\p n ) 
dqf s 



8=1 "Vis 8=1 

The force acts on the first particle. This force is 
a sum of the internal forces F^ k = F s (qi,pi,qfe,pfc,t), 
and the external force F} e = F s e (qi , px , t). In the case of 
binary interactions, we have 



^ = ^ le +E^ lfe - 



(91) 



fc=2 



Using Eq. I|91|) . the second term in Eq. I|9l)|l can be 

rewritten in the form 



10 



;« [2 ,...,„ l( ™ + g^) 



fc=2 



,«](^Pn). (92) 



Since the n-particle distribution function p n is a symmet- 
ric function for the identical particles, we have that all 
(n— 1) terms of sum l|92[) are identical. Therefore the sum 
can be replaced by one term with the multiplier (n — 1): 



S ;o ' 2 n] U F - 



P> 



(n-l)I a [2,...M ^(F^pn). 



Using 7 a [2,...,n] = I a [2]I a [3, ra], we have 



/ a [2]^(i^T[3,...,n]p n 



(9 



dp? 



-t[2] (f?~ P2 ). 



Here we use definition l|77|) of two-particle distribution 
function. Since pi is independent of q2, P2, we can 
change the order of the integrations and the differenti- 
ations: 



d P f s 



Finally, we obtain the equation for one-particle re- 
duced distribution function, 



dpi \ - v^lPD A ^ 
9t f^i dqf s + ^ 



8(0%) 



(n-l)AJ(p 2 ). 

(93) 



Here I{p~2) is a term with two-particle reduced distribu- 
tion function, 







2 P2- 



(94) 



Equation i|93|l is called a first fractional Bogoliubov equa- 
tion (first equation of Bogoliubov chain) . 

Let us consider the physical meaning of the term I(p2)- 
The term I(p2)dp a (q,p) is a velocity of particle num- 
ber change in 4m-dimcnsional elementary phase volume 
dfiaili, P2, q2, P2)- This change is caused by the inter- 
actions between particles. If a = 1, then we have the 
first Bogoliubov equation for non-Hamiltonian systems. 



C. Second Fractional Bogoliubov Equation 

The fractional Liouville equation allows us to derive 
equation for two-particle reduced distribution function 
p~2 in the form 



dp 2 
dt 



^A fej 5 2 + A 12 p 2 + c(n)AT( j 5 3 ), (95) 



where c(n) — (n — l)(n — 2)/2, and At is one-particle 
Liouville operator, 



d{G k s P2) 



A J2 



8=1 8=1 ^* 

and A12 is two-particle Liouville operator, 



m q m o 



and I(p~3) is a term with the three-particle, reduced dis- 
tribution 

s— 1 ^ s 



2 s 

The derivation of Eq. (|95|l is the analogous to the deriva- 
tion of Eq. 

It is easy to see that Eqs. I|93|) and (|95|l are not 
closed. The system of equations for the reduced dis- 
tribution functions are called the Bogoliubov hierarchy 
equations. 



VI. FRACTIONAL VLASOV EQUATION AND 
DEBYE RADIUS 



A. Fractional Vlasov Equation 

In this section, we derive the fractional analog of 
the Vlasov equation from the first fractional Bogoliubov 
equation. Let us consider the particles as statistical in- 
dependent systems. In this case, we have 

/5(qi,Pi,q2,P2,*) = |5(qi,Pi,t)p(q 2 ,P2,i)- (97) 
Substituting Eq. J§3 in Eq. l|9"I|). we get 

m ft 

r(P2) = -^^h[i]i a [2}Frp 1 [2]. 

Here we use the notation p[k] for the distribution function 
/o(qfe, pfe, t). As the result we have the effective forces, 

^ e// (qi,Pi,U = ^Wfpi[2]. 
In this case, we can rewrite the term 1)940 in the form 



8 



~(piF} eff ). 



I(P2 



Substituting Eq. (JHHJ in Eq. (EH), we get 



(98) 



fc=i 



dpi -A d{G\pi) d 
dt ^ dgf, ^ dp? 

S — l L * S—l ±lj 

(99) 

where b = n — 1. This equation is a closed equation for 
one-particle distribution function with the external force 
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F le and the effective force F le ^ . Equation l|99|l can be 
called a fractional Vlasov equation |37| . 

The fractional Liouville, Bogoliubov and Vlasov equa- 
tion are better approximation than its classical analogs 
for the systems with the fractional phase spaces (the frac- 
tal dimensional spaces). For example, the systems that 
live on some fractals (spaces with the fractal dimensions) 
can be described by the suggested fractional equations. 



B. Debye Radius 

In this section, we consider the Debye radius for the 
fractional systems that are defined by the equations 



ks 



pi 



dpi 



dt 



m 



dt 



(100) 



where we use the dimensionless variables qk, Pk, -Ffc, t. 
Let ro = qo be the radius of the interaction, Here m = 
Mr /toPo is a dimensionless mass, where M is a particle 
mass. Using m = 1, we get 



to 



Mqo 
Po 



A 



t F Mq Q 



Po 



Pi 



(101) 



Let us use the condition p\jM — kT$ for the character- 
istic momentum pg. Note that the condition p 2 /M = kT 
can be realized for non-Hamiltonian and dissipative sys- 
tems H3. 

The first fractional Bogoliubov equation 1)93(1 for the 
dimensionless one-particle distribution pi has the follow- 
ing dimensionless form |3S| : 



2£l + v d ^p^ 
at ^ 



dPis 



ABI(p 2 ). (102) 



The dimensionless first Bogoliubov equation l|102|) has 
two characteristic parameters, 



t F _ MrpFp _ roFo 



Po 



Pi 



B = n Q r { 



3<* 



(103) 



Let us consider the coefficient B. It is known that 
fractal particle system and fractal media are described by 
the power law relation (|29|l : 



nir) = nor 



D 



(104) 



where D < 3 and n is the D-dimensional concentration 
of the Z?-dimcnsional distribution of particles. The di- 
mension D of fractal system is connected with order of 
the fractional integrals a by D = 3a. The concentration 
no can be defined by the £>-dimensional mass density k: 
n = k/M which is used in Eqs. (|2Sl and (jSHJ- To cal- 
culate the mass fractal dimension D and concentration 
no, we can take the logarithm of both sides of Eq. (|104fl . 
When we graph ln(n) as a function of ln(r), we have 

ln(n) = D ln(r) + Zn(no), 



and we get a value of the fractal dimension D of fractal 
media and parameter n . Therefore these values can be 
measured for homogeneous fractal media. 

Let us consider the fractional systems with the force 



F, 



1 



ki 



47re ro \r k - n\ 



28 ■ 



(105) 



where and r; are dimensional values of coordinates. If 
5 = 1, then we have the usual electrostatic interaction. 
In this case, the Gauss theorem for the fractional space 
is not satisfied. If 2<5 = 3a — 1, then the Gauss theorem 
for the fractional space is satisfied. The radius ro and 
the force Fq are connected by the equation 



Fa 



1 



4tt£o 1q 



(106) 



Using the relation AB ~ 1, we have the characteristic 
radius of the interaction in the fractal media 



ro = r D 



_i e kT 



e 2 no 



(107) 



which can be called a fractional Debye radius. If the par- 
ticle systems or media have the integer mass dimension 
D = 3, then a = 1, and we get the usual equation for the 
Debye radius [3^ • The fractional radius characterizes the 
scale go = r o of the fractal media or fractal system with 
non-integer mass dimension. 



VII. CONCLUSION 

In this paper, we consider the fractional generalizations 
of the phase volume, the phase volume element and the 
Poisson brackets. These generalizations lead us to the 
fractional analog of the phase space. The space can be 
considered as a fractal dimensional space. We consider 
systems on the fractional phase space and the fractional 
analogs of the Hamilton equations. The physical inter- 
pretation of the fractional phase space is discussed. The 
fractional generalization of the average values is derived. 

It is known that the fractional derivative of a constant 
need not be zero. This relation leads to the correlation 
between coordinates q and momentum p. Therefore q and 
p are not independent variables in the usual sense. As 
the result, the generalization of the Poisson brackets with 
fractional derivatives l|15|) is not canonical. In order to 
derive equations with fractional derivatives we must have 
a generalization of Darboux theorem 40] for symplectic 
form with fractional exterior derivatives. However this 
generalization is an open question at this moment. In 
order to define Poisson brackets with the usual relations 
for the coordinates and the momenta we can use Poisson 
brackets (I22|l with the fractional power of coordinates 
and momenta. 

Note that the dissipative and non-Hamiltonian systems 
can have stationary states of the Hamiltonian systems 
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[4l|. Classical dissipative and non-Hamiltonian systems 
can have the canonical Gibbs distribution as a solution of 
the stationary Liouville equations for this dissipative sys- 
tem [5l| . Using the methods [2]| , it is easy to prove that 
some fractional dissipative systems can have fractional 
analog of the Gibbs distribution (non-Gaussian statistic) 
as a solution of the fractional Liouville equations. Us- 
ing the methods |2l|, it is easy to find the stationary 
solutions of the fractional Bogoliubov equations for the 
fractional systems. 

Note that the quantization of the fractional systems 
is a quantization of non-Hamiltonian dissipative sys- 
tems. Using the method, which is suggested in Refs. 
(42l E^L IH ] , we can realize the Weyl quantization for the 
fractional systems. The suggested fractional Hamilton 
and Liouville equations allow us to derive the fractional 
generalization for the quantum systems by methods sug- 
gested in Refs. HHSm. 

In this paper the fractional analogs of the Bogoliubov 



hierarchy equations are derived. In order to derive this 
analog we use the fractional Liouville equation we 
define the fractional average values and the fractional 
reduced distribution functions. The fractional analog of 
the Vlasov equation and the Debye radius are considered. 

The fractional Bogoliubov hierarchy equation can be 
used to derive the Enskog transport equation. The frac- 
tional analog of the hydrodynamics equations can be 
derived from the first fractional Bogoliubov equation. 
These equations will be considered in the next paper. 

It is known that the Fokker-Planck equation can be de- 
rived from the Bogoliubov hierarchy equations 8]. The 
fractional Fokker-Planck equation can be derived from 
the fractional Bogoliubov equation. However this frac- 
tional Fokker-Plank equation can be differed from the 
equation known in the literature 0, 0, ■ 

The quantum generalization of the suggested fractional 
Bogoliubov equation can be considered by the methods 
that are suggested in Refs. [43, LiJ ■ 
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